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Abstract: An analytical calculation of the interaction geometry of two 

interlinked second-order nonlinear processes fulfilling phase-matching 

conditions is presented. The method is developed for type-I uniaxial crystals 

^ . and gives the positions on a screen beyond the crystal of the entangled 

Si^ ' triplets generated by the interactions. The analytical results are compared 

^— ^ . to experiments realized in the macroscopic regime. Preliminary tests to 

tH- ' identify the triplets are also performed based on intensity correlations. 
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1. Introduction 

The production of multipartite entangled states by means of multiple nonlinear interactions oc- 
curring in a single nonlinear crystal has been recently suggested [1, 2, 3] as an alternative to the 
use of single-mode squeezed states [4, 5, 6] or of two-mode entangled states and linear optical 
elements [7]. As in the case of bipartite states produced by spontaneous parametric downcon- 
version, such multipartite states display entanglement also in the macroscopic regime in which 
bright outputs are generated: for this reason they are particularly interesting for all the applica- 
tions of continuous-variable entanglement [8, 9]. 

The possibility of realizing the simultaneous phase-matching (PM) of two traveling- wave para- 
metric processes in a single crystal in a seeded configuration has been already demonstrated 
[10, 1 1]. In this paper we present the experimental realization of the same interlinked interac- 
tions starting from vacuum fluctuations. The output of the nonlinear crystal displays the entire 
ensemble of inseparable tripartite entangled states satisfying the PM conditions. In order to 
identify a triplet, we develop a 3D calculation of the characteristics of interlinked PM that re- 
produce the experimental phenomenology. Finally we perform a preliminary verification of the 
number of photons conservation law by means of intensity correlations. 

2. Theory 

We consider the Hamiltonian describing the simultaneous PM of a downconversion and an 
upconversion processes 

^int = rifljfls + 720203+ h-C. , (1) 

where yi °^ 04 and 72 °^ 05 are coupling coefficients. The two interlinked interactions involve 
five fields Oj, two of which, say 04 and a^, enter the crystal and act as non-evolving pumps. 
When acting on the vacuum as the initial state, Hi„i admits the following conservation law 

Niit)=N2it)+N3{t), (2) 

being A^j(r) = («[ (f)flj(f)) the mean number of photons in the j-th mode, and yields a fully 
inseparable tripartite state [1]. We consider the realization of Eq. (1) in a negative uniaxial 
crystal in type-I non-collinear PM interaction geometry. The processes must satisfy energy- 
matching {CO4 = ffli + (B3, CO2 = CO3 + (O5) and PM conditions (^ = k" + kj, k, = k'^ + kj), 
where ©j are the angular frequencies, kj are the wavevectors and "■'' indicate ordinary and 
extraordinary field polarizations. In order to investigate the geometrical constraints imposed by 
the PM conditions, we analytically calculate the output angles of each interacting field. 

To simplify calculations, we assume, in the reference frame depicted in Fig. 1, that the two 
pumps lie in the {y,z)-plane containing the optical axis (OA) and the normal to the crystal 
entrance face and that the pump field 04 propagates along the normal, z- The solutions cor- 
responding to the effective experimental orientation of the crystal can be obtained by simply 
calculating the refraction of the beams at the crystal entrance/output faces. Accordingly, the 
PM conditions for the two interactions simultaneously phase-matched can be written as 

kisinfii+kismlii = (3) 

A:icos/3isiniJi +fc3Cosj33sini93 = (4) 

A:iCOSj8iCOSlJi +fc3COSj33COSj33 = ^4 (5) 

^2sin)32 = ^3sinj33 (6) 




Fig. 1. Scheme of the phase-matched interlinked interactions: (x,y)-plane coinciding with 
the crystal entrance face; a, tuning angle; jSj's, angles to (y,z)-plane; tJj's, angles on the 
(y,z)-plane; (p, angle to the optical axis (OA). 



ki cos j32 sin ■62 
^2 cos )32 cos 192 



^3 cos )33 sin dj+ks sin tJs 
^3 cos )33 cos ■d-T, + k^ cos d-^ 



(7) 
(8) 



where the angles are defined as in Fig. 1. The wave vectors k-^ are defined as k-^ = «j(a)j,kj)ci)j/c, 
being c the speed of light in the vacuum and «j(Oj,kj) the refraction indices of the medium 



Hj (cOj) = Ho (cOj) 



niioh,^) 



cos^^ 



sin^^ 
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-1/2 



;■ = 1,3,5 
n4(a)4,a) = 



cos^a 



sin^a 



«2(a)4) «2(a)4) 



-1/2 



(9) 

(10) 



where no,e{co) are given by the dispersion relations of the medium. Note that cos^ = 
cosj32 cos(t92 - «) and jis = ^4 = t}4 = 0. 

For fixed frequencies and propagation directions of the pump fields, we are left with 1 1 vari- 
ables (cOi, 0)2, (O3, tJi, '62^ t93, tJs, j3i, j32, Pi, (X) and 8 equations only (energy conservation and 
Eqs. (3) to (8)). The problem can thus be solved by choosing three of the variables (say (0\, TJ5 
and a) as parameters. The Appendix contains the algebraic procedure to analytically solve the 
problem. In order to efficiently handle the dependence of the solutions on the free parameters, 
we implemented the calculation by using the software Mathematica (Wolfram Research, IL). 

In Fig. 2 we show one picture of the movie containing the comparison between measured and 
calculated results for the output of the crystal as a function of the tuning angle a for a fixed value 
of the external angle between the two pumps. Note that not all the calculated output frequencies 
appear in the experimental part due to the sensitivity of the camera sensor. As the pictures of the 
experimental outputs were taken on a screen located beyond the crystal normally to the direction 
of pump fl4, the calculations had to take into account both the rotation of the crystal and the 
refraction of the beams at the entrance and exit faces of the crystal. In particular, refraction at 
the exit face gives: sinjSjout = njsinjSj and sinjJj out = «j/(l — «?sin2/3j)''2cosj3jSini?j. 

3. Experiment 

For the realization of the interaction described by Eq. ( 1 ), the pump fields were provided by the 
third-harmonics and the fundamental outputs of a continuous-wave mode-locked Nd: YLF laser 
regeneratively amplified at the repetition rate of 500 Hz (High Q Laser Production, Austria). In 




Fig. 2. (1.271 MB) Movie showing the variation of the measured and calculated outputs of 
the crystal as a function of the tuning angle a for fixed external angle between the pumps. 
Left panel: picture of the output of the crystal taken with a commercial digital photocamera 
(no color correction applied). LHS: a portion of the downconversion cones. RHS: output of 
the upconversion process. Right panel: calculated output of the upconversion process. 



particular, the third-harmonic field (X4 — 349 nm, ~ 4.45 ps pulse duration) was used as the 04 
field producing the downconversion cones and the fundamental field {X5 = 1047 nm, '^ 7.7 ps 
pulse duration) as the 05 field pumping the upconversion process. 
As depicted in Fig. 3 a), both pumps were focused and injected into a j3-BaB204 crystal (BBO, 
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Fig. 3. Upper panel: scheme of the experimental setup. BBO, nonlinear crystal; Fi_3, fil- 
ters; Pi_3, pin-holes; /1-5.S, lenses; Di_3, p-i-n photodiodes; SGI, synchronous gated- 
integrator; MCA+PC, computer integrated multichannel analyzer. Lower panel: picture of 
the visible portion of the output states on a screen located beyond the nonlinear crystal. 



Fujian Castech Crystals, China, 10 mm x 10 mm cross section, 4 mm thickness) cut for type- 
I interaction (jJcut = 34 deg) at the angle z?5_ext — —34.8 deg with respect to each other. The 
required superposition of the two pumps in time was obtained with a variable delay line. For 
alignment purposes, we used the light emitted by a He-Ne continuous-wave laser (Melles-Griot, 



CA, 5 mW max output power) to seed the process at Ai = 632.8 nm. The beam was colHmated 
and sent to the BBO in the plane containing the pumps at the external angle z?i,ext = —2.54 
deg with respect to 04. The seeded interactions produced two new fields: 03 (A3 — 11%. 2 nm, 
^3.ext — 3.35 deg) generated as the difference-frequency of a^ and a\, and 02 (■^2 = 446.4 
nm, i?2.ext = —12.78 deg) generated as the sum-frequency of 03 and 05. The pump fields were 
sufficiently intense so as to observe the process starting from vacuum fluctuations: as shown 
in the picture in Fig. 3 b), on a white screen located beyond the crystal it was possible to see 
both the tunable bright downconversion cones and the two polychromatic half-moons-shaped 
states generated by the upconversion process together with the spots of the fields generated by 
the seeded interactions. 

The conservation law in Eq. (2) implies strong intensity correlations among the generated 
fields which are necessary but not sufficient to demonstrate the entangled nature of the triplet. 
We realized intensity-correlation measurements by filtering each portion of light in frequency 
and aligning three pin-holes of suitable diameters on the spots of the seeded process. Note that 
the seeding He-Ne light was switched-off once completed the alignment. The light selected by 
the pin-holes was focused on three p-i-n photodiodes (two, Di 2 in Fig. 3 a), S5973-02 and one, 
D3, S3883, Hamamatsu, Japan). The overall detection efficiencies in the three detection arms 
were Tji = 0.44, 772 = 0.72 and T73 = 0.43. Each current output was integrated by a synchronous 
gated-integrator (SGI in Fig. 3 a)) in external trigger modality, digitized by a 13-bit converter 
(SR250, Stanford Research Systems, CA) and recorded in a PC -based multichannel analyzer 
Data acquisition and analysis were performed with software Lab View (National Instruments, 
TX). We evaluated the following correlation function on subsequent laser shots (k) 

^ ^ ([OTi(i) - (mi)] [OT2(i + k) +OT3(i + k) - {1112 + 1113)]) 
C7{mi)a{m2 + mj,) 

in which otj is the number of detected photons and a{x) = ((x^) — (x)^)''^ is the standard 
deviation. It can be demonstrated [12] that in the case of bright triplets the correlation coefficient 
e = r(0) must approach unity. The mean values of the detected photons in the three parties of 
the triplet were Mi = 1.082 x 10^ M2 = 2.3 x 10'' and M3 = 1.115 x 10^ and the measured 
correlation coefficient was e = 0.916. 

4. Conclusions 

The perfect agreement between calculations and experimental results demonstrates the correct- 
ness of the developed model. This enables us to forecast different interaction schemes that in 
the future will allow to overcome the experimental difficulties experienced in the present setup. 
In fact, the reason why the measured correlation coefficient is less than unity could be the 
imperfect selection of the triplet and/or the presence of spurious light in the same location. 
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Appendix 

We describe the analytical procedure to solve the equation system (3)-(8). 



From the first of Eqs. (10) we write 

A:2= [G2Cos2(j32-a)cos2j32+L2]"'/^ , (12) 

where G2 = c^/co^[l/}^l{o}2) -l/njiah)] = 1/^2,0^ lAie and L2 = c^/(oll/nl{o}2) = 
1/^2 g. By squaring and summing Eqs. (3), (4) and (5) and defining A = {k^ + k^ — k^)/{2kik4), 
we easily get 

COSJ33-— V' (13) 

COSO3 

whereas by squaring and summing Eqs. (6), (7) and (8) and Eqs. (7) and (8) we get 

kj = A:3 +A;5+ 2/^3/^5 cos )33 (sin )33 sin i?5+ cos ft cos iJs) (14) 

/fejCOS^ft = A:3COS^j33+/t5+2/t3/t5COS)33(sin)33SiniJ5+COSj33COStJ5) . (15) 

By eliminating ^2 from Eqs. (14) and (15) and using Eq. (13) we find 

, „ klA^ tan^ ^^3 + yt? + Ik^k^A (tan ^^3 sin t^s + cos tJs ) 

cos JSt = ^ T T ^ ^ ^ . (16) 

/tf + ytf + 2/t3yt5A (tan tJ3 sin 1^5 + cos 1^5 ) 



From Eq. (3) we obtain 



sin 61 = --—sin 133 . (17) 

ki 



that, once substituted into Eq. (4) together with Eq. (13), gives 



k^A 
sintJi = , (18) 



//t?-/t2sin2j33 

We now divide Eq. (7) by Eq. (8) and use Eq. (13) 

^3Atant93+fc5sintJ5 

tani32 = — , ■ , , ^^ • (19) 

K3A + K5 cos v-^ 

By squaring Eq. (12) and inserting Eq. (14) and Eq. (16) into it we get 

2 1 - L2 [^3 +kl + Ikjks cos fe (sin fe sin ^s + cos fe cos i^s)] 

G2 [k^ cos^ /33 + ^5 + 2^:3^5 cos /33 (sin j33 sin d-s + cos ft cos 7^5 )] 

By substituting Eq. (13) into Eq. (20), expanding cos^ {d2 — oc) in Eq. (20) and using Eq. (19) 
to eliminate i!h, after some simplifications we get a quadratic algebraic equation for the variable 
tan 1^3 

atan^j33+^tanjJ3+c = 0, (21) 

whose coefficients are 

a = G2A ^3 sin' a 

b = 2G2A kT, (A kT, + k^ cos d^ ) sin a cos a + 2G2A ^3^5 sin d-^ sin' a + 2L2^ ^3^5 sin d-s 

c ~ G2(A ^3+^5CosiJ5)^(cos^a-sin^a)+2G2fc5siniJ5(A ^3 + ^5 cos zJs) sin a cos a 

+ G2(A^/t3+/t5+2Ayt3yt5COSTJ5)sin'a+L2(/t3+yt5+2Ayt3yt5CoszJ5)-l . 

Once solved Eq. (21), we find d^ and then all the other variables as a function of the free 
parameters coi, a and 1^5. 



